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Abstract. We review mathematical results concerning exponentially small 
corrections to adiabatic approximations and Born-Oppcnhcimer approxima- 
tions. 



Introduction 

The goal of this paper is to review recent results on exponentially small non- 
adiabatic transitions in quantum mechanics. In Section 1, we provide background 
information about adiabatic approximations. In Section 2, we discuss the determi- 
nation of non-adiabatic scattering transition amplitudes. In Section 3, we describe 
the time development of exponentially small non-adiabatic transitions. In Section 
4, we turn to exponentially accurate Born-Oppenheimer approximations. Then 
in Section 5, we discuss the determination of non-adiabatic corrections to Born- 
Oppenheimer approximations in a scattering situation. 

1. Adiabatic Background 

The adiabatic approximation in quantum mechanics concerns the time depen- 
dent Schrodinger equation when the Hamiltonian depends on time, but varies on a 
very long time scale. Mathematically, this situation corresponds to the singularly 
perturbed initial value problem 

(1.1) tedtMt) = H(t)ip e (t), V 6 (0) = 0(0), 

where t £ R, e is a small parameter, and ip e (t) belongs to a separable Hilbert space 

n. 

In a simple situation, the adiabatic approximation relies on two basic assump- 
tions. The first is a regularity condition on the Hamiltonian: 

R: The Hamiltonian H(t) is a bounded self-adjoint operator on Ti that depends 
smoothly on t. 
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We denote the unitary propagator corresponding to (|1 . If) by U e (t, s). It satisfies 

(1.2) ied t U e (t )S ) = H(t)U e (t,s), U € (s,s) = I 
for all (t, s) £ R 2 

The second assumption is that there is a gap in the spectrum of H(t) at all 
times t: 

G: The spectrum of H(t), denoted by a(t) C R, can be decomposed as 
cr(t) — oi{t) U o-2(t), with inf dist(ai(t), <72(t)) = g > 0. 

Under this hypothesis there exists a self-adjoint spectral projection P(t) corre- 
sponding to o~ i (i) , which depends smoothly on t. 

In this situation, the Adiabatic Theorem of Quantum Mechanics says the fol- 
lowing: 

Theorem 1.1. Let H(t) satisfy R and G, and let P(t) be as above. Then, for 
small e, 

(1.3) \\(l-P(t))U e (t,s)P(s)\\ = 0(e), 
for any (t,s), where the error term is sharp. 

In other words, for small e, the quantum evolution follows the isolated spectral 
subspaces of the Hamiltonian, up to an error of order e. We call the error term the 
transition amplitude between the subspaces P(s)Tt and (I — P(t))H.. 

One modern proof actually yields a stronger result. Let us sketch it for later 
reference. Let V(t, s) be the (e-dependent) unitary evolution defined by 

(1.4) ied t V(t,s) = {H{t) - ie[d t P{t), P(t)])V(t,s) and V{s,s) = I. 
It is not difficult to verify the intertwining property 

(1.5) V(t,s)P(s) = P(t)V(t,s), for all (t, s) e R 2 . 

A generalized integration by parts procedure then establishes the estimate || U e (t, s) — 
V(t, s)|| = 0(e\t — s\), which implies the result. 

Historically, the Adiabatic Theorem goes back to Born and Fock j£j who proved 
it in 1928 for matrix valued Hamiltonians with non-degenerate eigenvalues. In the 
fifties, Kato [31] proved the statement for general Hamiltonians with an isolated 
eigenvalue E(t). Thirty years later, Nenciu showed the isolated eigenvalue E(t) 
could be replaced by an isolated piece of spectrum o~\{t). The extension to un- 
bounded operators was proven by Avron, Seiler and Yaffe in |5]. 

Although we shall always assume the gap hypothesis G, there are recent ver- 
sions of the adiabatic theorem that hold without it. Instead, these results require 
the existence of a sufficiently regular, finite rank family of projectors P(t) that sat- 
isfy P(t)H(t) = E(t)P(t) for all t. Here E(t) is an eigenvalue of H(t), and P(t) is a 
spectral projector of H(t) for almost all t. See e.g., |10| . [J, |43| . In this situation, 
the transition amplitude can go to zero arbitrarily slowly as e tends to zero. 

Since the early thirties, it has been known that the adiabatic theorem could 
be substantially improved under certain circumstances. Zener |45l considered p. 1(1 
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for the 2x2 matrix Hamiltonian 

(1-6) H{t) = \ (j ! t ), 

whose spectrum is a (t) = {±\ Vt 2 +S 2 }. For this Hamiltonian. fTTT|> can be 
solved exactly in terms of parabolic cylinder functions, for any 6 and any e > 0. In 
particular, in the scattering regime where the initial and final times s and t tend 
to — oo and +00 respectively, the transition amplitude takes the simple form 

(1.7) A(e) := ^lim^ || (I — P(t)) U e (t, s) P(s) || = e ^^l {At \ 

Zener's analysis shows that the transition amplitude decreases from its order e 
value for finite times (s,t) to an exponentially small value in 1/e in the scattering 
limit. Shortly after Zener's result, Landau |32j argued that l|1.7|l should also hold 
for more general analytic 2x2 real symmetric Hamiltonians whose non-degenerate 
eigenvalues Ei(t) and E 2 (t) displayed an avoided crossing with gap S, i.e., when 

(1.8) E 2 (t)-E 1 (t) ~ y/(t-t ) 2 + S 2 . 

Here 5 is a small, but fixed, extra parameter in the problem. These considerations 
gave rise to the famous Landau Zcncr formula, which showed that non-adiabatic 
transitions between spectrally isolated subspaces belonged to the realm of expo- 
nential asymptotics. 

This Landau-Zener mechanism is the main concern of this review. 

2. Exponentially Small Transitions in the Adiabatic Approximation 

From a mathematical point of view, the first exponential bounds on transition 
amplitudes in a general framework were obtained in |27j . If we assume for simplic- 
ity that the Hamiltonian is bounded, this result requires the following regularity 
and scattering hypotheses: 

A: There exists fx > such that the map 1 1— > H(t) £ £(7~t) has an analytic exten- 
sion to the strip = {z £ C : |/m(z)| < /x}. 

S: There exists v > 1, two bounded self-adjoint operators H(±oo), and a constant 
c, such that as t — * ±00, 

sup \\H(t + i8)-H(±oo)\\ < 

\s\<n W 

where (t) = yl + t 2 . 

Theorem 2.1. [27] Assume H(t) satisfies A, S and G. Then, there exist C 
and r > such that A(e) < Ce~ r / C , for small e > 0. 

The original proof in |27| establishes properties of solutions to the Schrodinger 
equation for complex t, by using complex WKB techniques. Similar results were 
subsequently proven using other methods. See e.g., |38) . |41| . |28j . and |35| . 

A rigorous asymptotic formula for the transition amplitude across the gap of 
a two- level system was proven in [25] . In addition to the assumptions above, the 
proof required supplementary hypotheses on the behavior in the complex plane of 
the so-called Stokes lines of the problem. 
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This hypothesis, that we call SL, is complicated, and we shall not attempt to 
describe it here. It is typical of the complex WKB method and is discussed in detail 
in [23] . 

Theorem 2.2. Let 



be a 2 x 2 hermitian matrix that satisfies assumptions A, S and G. Further as- 
sume that the analytic extension of p(t) = ^JX 2 {t) + Y 2 (t) + Z 2 (t) has only one 
conjugate pair of zeros, {zq, zq}, in and that the supplementary hypothesis SL 
is satisfied. Then, there exist G ^ and 7 > 0, such that, as e — > 0, 



Remarks 2.3. 1) The point zo is a complex crossing point for the analytic 
continuations of the eigenvalues of H(t), and the decay rate 7 is the one predicted 
by Landau. 

2) The prefactor G is equal to one for generic real- valued symmetric Hamiltonians. 
It can take other values for generic hermitian matrices, as was independently dis- 
covered by Berry |3] . 

3) Several variants and improvements of this theorem have been proven. See |22| . 
|28| . and p3 . Extensions to more general scattering systems were proved in 29 , 
|34| . and |24| . In particular, |23l I24| show that in an avoided crossing regime, 
hypothesis SL is automatically satisfied. See also |30| . 

These results were first proven by means of the complex WKB technique. An- 
other popular and fruitful method in exponential asymptotics is the so-called opti- 
mal truncation method. In this method, one first derives an asymptotic expansion 
of the desired quantity in powers of e. Next, one carefully estimates the difference 
between the exact quantity and the truncation of the expansion after n terms. This 
error estimate depends on both n and e, and if it is bounded by cq c™ n\ e™, one can 
choose n to depend on e in an optimal way to obtain an error estimate that is 
0(e~ v ( g '< e ). The optimal n satisfies n ~ g/e, where g > is sufficiently small. 

Nenciu 38 first implemented these ideas in a general adiabatic context. His 
work followed contributions by Garrido (12) and Sancho (40) in the sixties. The 
optimal truncation technique was also used in |28j , and an elementary exposition of 
it is provided in |19j . More recently, |39| and |44] further adapted and generalized 
these ideas to deal with space-time adiabatic theory. 

Let us briefly present this circle of ideas through the construction used in 28 . 
For a Hamiltonian Hit) satisfying R and G, we set H (t) — H(t), Po(t) = P(t) and 
Ko(t) = [dtPo(t), Po(t)]. Then, we iteratively construct a sequence of self-adjoint 
operators as follows. For any integer q > and e sufficiently small, we define 




X(t) + iY(t) 
-Z(t) 





(2.2) 



H q (t,e) 
K q (t,e) 



H(t) - ieK q ^(t,e), 
[dtP q (t,e), P g (t,e)], 



and 



where P q (t, e) is the spectral projector of H q (t, e) corresponding to Pq in the limit 
e — > 0. Perturbation techniques ensure that this scheme is well-defined for any q, 
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provided e is small. Furthermore, hypothesis S implies 
(2.3) Pq(t,e) = P a (t) + 0(e), as e -> 0, and 

lim P„(t,e) — P(±oo), for any fixed e. 

We next introduce the evolution V^(£,s) associated with H q (t,e) as the solution 
to <|1.4|) . with iJ g in place of -ff and K q in place of -Ko- By construction, the 
intertwining relation 1)1. 5|1 holds with V q and P q in place of and P. The motivation 
for this construction is that it yields the estimate 

\\U e (t,s)-V q (t,s)\\ < J \\K q (u,e)-K q -i(u,e)\\du < \t - s\ f3 q e q+1 , 

for small e, where (3 q is finite. When H{t) satisfies the analyticity hypothesis A, 
one can further control (3 q as a function of q and prove the bound /3 q < Co c\ q\. The 
optimization over q can then be performed. This implies the following theorem. 

Theorem 2.4. [3H Assume H(t) satisfies A and G. There exist T > 
and on orthogonal projector P*(i, e) = P(i) + 0(e) suc/i t/iat 

||(I-P*(i,e))C/ £ (*,s)P*(s,e)|| = 0(e- r / £ |i-s|). 

Remarks 2.5. 1. This result shows that although the transitions between 
instantaneous eigenprojectors of H(t) are of order e for finite times, there exists 
another set of orthogonal projectors which are better suited to the adiabatic ap- 
proximation, in the sense that transitions between them are always exponentially 
small. 

2. We call the projectors P*(i, e) optimal adiabatic projectors. They are not 
uniquely defined. 

3. We recover Theorem 12. II if we further assume S. Sec 1)2. 3[l . 

4. The construction above was used in |28) and |23) to reduce the study of tran- 
sitions between two isolated levels in the spectrum of a general Hamiltonian to an 
effective reduced 2-level problem. 

3. Time Development of Exponentially Small Non Adiabatic 

Transitions 

Berry [5] pushes these ideas further for two-level systems. For such systems, 
a one-dimensional projector corresponds to an orthonormal basis of C 2 , with one 
vector in the range and one vector in the kernel of the projector. Berry proposes 
studying the time development of the transition amplitude in this time-dependent 
basis. His heuristic arguments [2] are supported by spectacular numerical results 
|33| . These show that the transition amplitude passes from its initial value of zero 
to its final exponentially small value l|2.1|l in a universal way described by an error 
function. 

The first rigorous results 20 on this question were obtained for the particular 
2x2 Hamiltonian function 

( "» "«> " uWTW ( < ) ■ 

for fixed S > 0. The eigenvalues are E\(t) = 1/2 and E^t) = — 1/2 for every t. 

This Hamiltonian can be viewed as the Landau-Zcncr Hamiltonian modified 
to keep its eigenvalues constant and with the gap normalized to the value 1. 
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The notion of avoided crossing has been replaced by the singularities of the 
Hamiltonian at t = ±iS. These points govern the transitions between the two 
levels and yield the transition amplitude y/2e~ s ^ e . See |22| . The time development 
is described by the following result. 

Theorem 3.1. |20| Let H(t) be given bv $3.1]) and let u>i (t) andip 2 (t) be smooth 
real-valued normalized eigenvectors corresponding to the eigenvalues 1/2 and —1/2, 
respectively. Then for any fi S (0, 1/2), 

1. There exist vectors Xi(e, t) and X 2(e, t) that satisfy the Schrodinger equation 
11. 1)) up to errors of order e -<5//e and correspond to the eigenstates <fi\(t) and (f2(i) 
in the sense that 

(3.2) Jim | (<pj(t), Xi (e, t))\ = 1 + 0{e~ s '^). 

|t|— >oo 

Moreover, the set {Xj( e i t)}j=i,2 is orthonormal up to errors of order e~ s ^ . 

2. The Schrodinger equation has solutions ^j(e, t), j = 1,2, such that uniformly 
in t G M as e — > 0, 

*x(e, t) = xi(e, t) - V2e- 3 /* 1 1 erf + l} Xa(e, *) + 0( e - 5 / e e ^), 

and 

* 2 (e, t) = X2 (e, t) + V2 e- s ^ \ jerf (^=) + 1 } Xi{e, t) + O^e^). 
Remarks 3.2. 1. Recall that the function erf is defined by 

(3.3) erffx) = -= / e~ v% dy e [-1, 1]. 

V 71 " Jo 

The vectors Xj( e i J = 1> 2, are explicitly constructed as approximate solutions to 
(11.11) obtained by means of optimal truncation of asymptotic expansions of actual 
solutions. As t — > — oo they are asymptotic to the instantaneous eigenvectors ifj(t) 
of H(t), up to a phase. We call them the optimal adiabatic states. 

2. The transition mechanism between optimal adiabatic states goes from the value 
zero to the value y/2 e~ s / e in a smooth monotonic way described by the switching 
function (erf+l)/2, on a time scale of order -Je. By contrast, the transition between 
instantaneous eigenstates of the Hamiltonian displays oscillations whose amplitudes 
are of order e for any finite time. It reaches its exponentially small value only at 
t = oo. The vectors Xj'( e i t) are a ls° optimal in that sense. 

3. Since the optimal adiabatic states and eigenstates essentially coincide at t = 

±oo, the transition amplitude equals -4(e) — \/2e~ s ^ e , up to errors of order 

e -5/e £ ^ 

4. The whole evolution operator U e (t, s) associated with 111 .21 is actually known 
up to errors of order e _<5 / e e**, for any t and s. 

We prove this theorem by carefully controlling the terms that occur in the a- 
symptotic expansion of exact solutions, so that we obtain very detailed information 
about the optimally truncated expansion. Even for our simple Hamiltonian, this in- 
volves very delicate asymptotic analyses of solutions to several non-linear recursion 
relations. 

This result was soon followed by an essentially equivalent result by Betz and 
Teufel 5 , who dealt with a slightly larger family of Hamiltonians. The main 
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point of their analysis was to introduce a different technique. They defined the 
optimal adiabatic states as instantaneous eigenstates of transformed Hamiltonians, 
similar to l|2.2|) . instead of dealing with the expansions of solutions. This way they 
simplified the analysis of the recursion relations. That paved the way for their 
significant generalization 6 . They proved that the conclusions of Theorem 13.11 
held for general 2x2 time dependent analytic Hamiltonians. 

Let us informally state their result for real symmetric matrices 

Z{s) X{s) \ 
X(s) -Z(s) ) ' 

such that p(s) = y / X 2 (s) + Z 2 (s) is strictly positive for real s and has only one 
zero zo G C + , the upper half plane. Under certain technical assumptions, Betz and 
Teufel show that the conclusion to Theorem 13. II holds with t given by 

t( s ) — / p{u)du, 
Jo 

5 given by 5 — |Im t(s)\ s=zo \, and the prefactor v2 in front of the curly brackets 
replaced by G ^ 0, which depends on the properties of H (s) in a neighborhood of 

Z . 

See [5] for a non-technical account of their findings, which confirm Berry's 
intuitive arguments. 

4. Exponentially Accurate Born Oppenheimer Approximation 

We now turn to transition phenomena in the Born-Oppenheimer approxima- 
tion. This approximation is widely used in molecular dynamics. It takes advantage 
of the disparity between the masses of nuclei and electrons, and accordingly, of the 
very different time scales for these particles' motions. The light electrons adjust 
quickly to the relatively slow changes in the configuration of the heavy nuclei. If 
they start in an electronic bound state, they stay in that bound state, as if the 
nuclei were not moving. This is typical adiabatic behavior. The nuclei are heavy, 
so they behave in a semiclassical fashion. The two motions are coupled because the 
electronic energy level plays the role of an effective potential for the nuclei. 

Hence, the Born-Oppenheimer approximation is simultaneously adiabatic and 
semiclassical in nature. The breakdown of these approximations are both associated 
with exponentially small phenomena: non-adiabatic transitions on the one hand, 
and tunnelling on the other hand. See e.g., |16L 1171 111) . It is thus reasonable to 
expect the Born-Oppenheimer approximation to be valid up to exponentially small 
corrections. 

We set the electron masses equal to one and take the nuclear masses propor- 
tional to e~ 4 , where e is a small parameter. We let x E M. d denote the collective 
positions of the nuclei and y £ l n denote those of the electrons. If W(x,y) is the 
inter-particle potential, the molecular Hamiltonian is 

(4-1) H{e) = -^-A x -±A y + W(x,y), 

on the Hilbert space L 2 (M. d x R"), where A z denotes the Laplacian in the z vari- 
ables. We introduce the electronic Hilbert space Ti e i = L 2 (W l ) and the electronic 
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Hamiltonian h{x) — — \ A y + W(x, y) on TC e i, which depends parametrically on 
the position x of the nuclei. Then we can rewrite l|4.1|l as 

(4.2) H{e) = -- A x + h(x) 

In an appropriate time scale, the molecular time-dependent Schrodinger equation 
is 

A 

(4.3) ie 2 d t V = -yA x $ + h{x)V, * £ L 2 (R d ,H e i). 

We study solutions to this equation in the small e limit, for times t £ [0, T], where 
T is fixed. We never actually use the explicit form of the Hamiltonian 1)4.1(1 . We 
study solutions to 1(4.3(1 with 14.2(1 satisfying the following two hypotheses. The 
first one says x <— > h(x) is analytic in an appropriate sense: 

A e ;: (i) For any x £ M. d , h(x) is a self-adjoint operator on some dense domain 
T> C TLeh where TL e i is the electronic Hilbert space. We assume the domain T> is 
independent of x, and that h(x) is bounded below, uniformly for x £ M. d . 

(ii) There exists fi > 0, such that for every ip £ T> ; the vector h(x)ij) is analytic 
in = {z £ C d : \Im(zj)\ < (i, j = l,...,d}. 

The second hypothesis is the familiar gap hypothesis, expressed this time for a 
space variable: 

G e i m . There exists an open set 5 C R d , such that for all x £ S 7 there exists 
an isolated, multiplicity one eigenvalue E{x) of h(x) associated with a normalized 
eigenvector $(x) £ Ti. e i- 

Under these two hypotheses, the Born-Oppenheimer approximation colloqui- 
ally says the following: Pick an initial condition for 1(4. 3() that is the product of a 
semiclassical wave packet concentrated near (a(0),r/(0)) in phase space, times the 
electronic eigenstate $(x). Then, at some later time t > 0, the solution \1/ is again 
given to leading order in e by a semiclassical wave packet times the same electronic 
eigenstate. The nuclear wave packet is now concentrated near the phase space 
point (a(t), T)(t)), determined from (a(0),ry(0)) by the classical flow generated by 
the classical Hamiltonian p 2 /2 + E(x). The spreading of the nuclear wave packet 
is also determined by the classical dynamics. 

As implied above, there is an exponentially accurate version of the Born- 
Oppenheimer approximation. 



Theorem 4.1. |18) Assume hypotheses A e i and G e ;. There exists a vector 
**(a;,t,e) £ L 2 (R d , TL e i) (that depends on a parameter g) for t £ [0, T], such that 
for small values of g, there exist C(g) and T(g) > 0, such that for small e > 0, 

e -« B W/ 1> $,(i 1 ) £) - *,(x,t,e) < C{g) e~ T ^/ e2 . 

In the state ^ , *(x,t,e), the electrons have a probability of order e 4 to be in the 
electronic subspace Span{^>(x)}- L . For any b > and sufficiently small values of g, 
the nuclei are localized near a classical path a(t) in the sense that there exist c(g) 
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and 7(g) > 0, such that for small e > 0, 



The nuclear configuration a(t) is determined by classical dynamics in the effective 
potential E{x). The limiting time T > is determined by the condition {a{t), t € 



Remarks 4.2. 1. Under more restrictive hypotheses, the limiting time T can 
be allowed to grow like | ln(e 2 )|, the Ehrenfest time scale. The error term becomes 
e -r(a)/£" ) where o < cr < 2. 

2. The exponentially accurate approximation has non-trivial components in the 
electronic space orthogonal to the electronic eigenvector <fr(x) which are of order e 2 , 
but its semiclassical dynamics is entirely determined by the electronic level E{x). 
These components of order e 2 correspond to the adjustment of the spectral subspace 
necessary to achieve exponential accuracy at finite times. See Ij2.3|) . 

We shall not be more precise here about the construction of the approximate 
solution We do not need it, and the details can be found in |18j . We just men- 
tion that $ * is obtained by optimal truncation of an asymptotic series in powers 
of e for the solution to l|4.3() . Results of the same sort have then been obtained in 
|39L 1361 142 j in more general situations and in other contexts. Roughly speaking, 
these papers show that transitions between (adjusted) spectrally isolated electronic 
subspaces are exponentially small under the time evolution. They can handle nu- 
clear wave packets which are not necessarily concentrated near a point of phase 
space in the semiclassical limit. In order to determine the characteristics of the 
nuclear component of the solutions, a semiclassical analysis within the spectral 
subspaces needs to be done. An important feature of these papers is that they 
decouple the adiabatic and the semiclassical effects. 

In any case, we infer from these results that transitions between electronic 
subspaces are hidden in the exponentially small error term. They are expected to 
stem from some Landau-Zener mechanism. The results in the next section verify 
that this expectation is correct. 

5. Non- Adiabatic Corrections to the Born Oppenheimer 

Approximation 

In Theorem 14.11 the nuclear wave packet is concentrated near a(t). Thus, it is 
reasonable to expect the electrons to behave as though their motion were governed 
by h(a(t)). Furthermore, it is reasonable to expect the Landau-Zener mechanism 
to govern their transitions, since the replacement of h(x) by h(a(t)) leads to a 
standard adiabatic approximation. 

From the analysis |21| of a simple Born-Oppenheimer model that we discuss 
below, the main qualitative features of this argument are correct. However, the 
quantitative predictions based on the use of the Landau-Zener formula are wrong 
at leading order. For example, the exponential decay rate is wrong, even if the width 
of the nuclear wave packet is negligible in the limit e — > 0. Some characteristics of 
the energy/momentum density of the nuclear wave packet are present in the main 
features of the piece of the wave function that has made a transition. These aspects 
are clearly absent from the heuristic considerations above. 





10 



GEORGE A. HAGEDORN AND ALAIN JOYE 



The following situation is considered in j21j : The nuclear variable x belongs 
to K and the electronic Hamiltonian is given by a finite order matrix h(x,d). We 
restrict attention to a 2 x 2 matrix for simplicity. We assume this matrix depends 
analytically on x and approaches limits sufficiently rapidly as x — > ±oo. As we 
saw previously, this is the most convenient framework in which a Landau-Zener 
formula can be proven. We finally assume that the two non-degenerate electronic 
levels display only one avoided crossing in a neighborhood of (x, 6) = (0, 0), of the 
type described by l|1.8f) . Here again, 6 is a supplementary parameter of the prob- 
lem. More precisely, with the notation of A e i for d — 1, we assume: 

M e ;: i) For some Sq > and /i > 0, the map x [0, Sq] 9 (x, S) i— > h(x, S) € M2(C) 
is C 2 . For any fixed S, x > /i(a;, (5) is analytic in S^. 

ii) There exist v > 5/2, c > 0, and iwo matrices h(±oo,5), such that 

sup ||/i(a; + iy, 5) — /i(±oo, (5)11 < , as ir — > ±oo. 

li/l<M.*e[o,<s ] \ x / 

The limiting matrices /i(±oo, 5) are C 2 in S G [0, <5o] and are non-degenerate. 

For each j£l and each 5 € [0, oo], we denote the two eigenvalues of h(x, S) by 
E\(x,S) and E2(x,S). We assume: 

AC e ;: i) For S > and x € [— oo, +oo], f/ie eigenvalues satisfy Ei(x,6) < E2(x,5). 

ii) When 5 = 0, the analytic eigenvalues Ej(x,0) have only one real crossing 
at x — 0. TTiey are labeled so that Ei(x,0) < E2(x,0) for all x < 0. 

iii) Moreover, we assume the following behavior near the avoided crossing: 

E 2 {x 1 5)~E 1 (x,S) = ^x 2 + S 2 +R 3 {x,5), 
where R3 is a remainder of order 3 around (0, 0) 

The corresponding normalized eigenvectors are denoted by <f>\(x,$) and 4>2(x,5). 
They are chosen to satisfy the phase condition (<f>j(x,d), d x <fij(x,8)) = 0. 

Remark 5.1. In the following discussion, 6 is positive and fixed, so we often 
drop it from the notation. We reintroduce it when necessary. 

The analysis of 1)4.3(1 in this set-up starts with a separation of variables and the 
study of generalized eigenvectors <!> of H(e). These vectors depend on a real energy 
parameter E, and they are defined by the ODE in C 2 : 

(5.1) (-J d * + $ = E *- 

We pick E in some bounded energy interval Act. To ensure that scattering is 
possible, we assume 

(5.2) inf A > sup E 2 {x). 

EeA xm 

Choosing an energy density Q(E,e) € C, we get an exact solution to l|4.3|l of the 
form 

(5.3) / Q{E,e)<$>(x,E,e)e- ltE /* 2 dE. 

J A 
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From physical intuition, we anticipate that if the energy density is sharply 
peaked, then when \t\ is large, these solutions to l|4.3|) will be concentrated in a 
region where \x\ is large. 

We choose the energy density to behave like a Gaussian of width e 2 , centered 
around an initial energy Eq in the interior of A. More precisely we assume: 

ED: The complex valued energy density is supported on A and has the form 
Q(E, e) = e- G ^ e2 e~ iJ ^ P(E,e), where 

i) G S C 3 (A) is non-negative, and satisfies G(E) = g(E-E ) 2 /2 + O(E-E ) 3 , 
for some Eq in the interior of A and some g > 0. 

ii) J G C 3 (A) is real-valued. 

iii) P(E,e)eC 1 (A) satisfies sup BeAje > \d%P(E,e)\ < C n , for n = 0,1. 

This choice implies that in the remote past or remote future the nuclear com- 
ponents are coherent states which basically propagate freely along the different 
electronic levels. Modulo the technicality of inserting a cut-off, Gaussian coherent 
states have such an energy density. For these states we use the notation of |15| : 

1 / B (x — a) 2 7] (x — a) \ 

(5.4) <po(A,B,e ,a,r,,x) = exp^- + i - , 

where the complex numbers A and B satisfy the normalization condition Re BA = 
1. These states are localized in position near x = a, and in momentum near p = r\. 
Their position uncertainty is e\A\ and their momentum uncertainty is e\B\. For 
a thorough discussion of these wave packets, see |15| . Note, however, that there 
are states characterized by a Gaussian energy density which are not necessarily 
Gaussian in space and momentum. See Section 6 of |21| . 

The solutions of 1|4.3[) considered in |21| are characterized by their behavior 
in the remote past. We now describe these asymptotic states. The generalized 
eigenvector equation H5.1(l can be solved by means of a WKB Ansatz as follows: 
Let kj(x,E) = \J2(E — Ej(x)), j = 1,2, be the classical momenta corresponding 
to the potentials Ej(x). Then, the solutions &(x,E, e) to the stationary equation 
(|5.1(l can be written as 

(5.5) *(x,E,e) = £ M*) 4ftE4 e-Wo-^W", 

3=1,2 \/£Kj{x,h) 



where Cj(x,E,e) € C are coefficients that satisfy some linear ODE in x, for each 
E e A. The condition M e ; ensures the existence of the limits of <pj (x) , kj (x,E), and 
dj(x,E,e), asu ±oo and further provides rates at which they are approached. 

We assume without loss that the "initial" conditions for the coefficients cj are 
fixed at x — ™oo and that cj(— oo, E, e) are independent of e and E. It follows that 
the dependence in E of all the above quantities is analytic in a neighborhood of A, 
and dj(x, E, e) is uniformly bounded in (x, E, e). 
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We now introduce freely propagating states ^"X^i £j e ; ±) E L 2 (M, C 2 ) that 
describe the asymptotics of the solutions ^(x, t, e) of (|4.3|l as t — > ±00. Let 

(5.6) i>°{x,t,e,±) = E rl>j(x,t,€,±) 

.3=1,2 

,^V JV ; i A A /2fe,-(±oo,£;) 



.3=1,2 



where cj^ (i?) = lim^^oo Jq 2 kj (y, E) dy~xkj (±00, E). These states are linear com- 
binations of products of free scalar wave packets in constant scalar potentials times 
eigenvectors of the electronic Hamiltonian. Their propagation is thus governed by 
the various channel Hamiltonians. With our sign convention, the states indexed by 
a = — propagate to the right, whereas those indexed by a — + travel to the left. 

By construction, ip~ (x, t, e, ±) is likely to be a good approximation of 'J (2, t, e) 
only for x — > ±00, i.e., for t — > ±00. Similarly, ip + (x, t, e, ±) should be a good 
approximation only for x — > ±00, i.e. for t — > =Foo. Quantitatively, we have 

Proposition 5.2. Assume AC e! , M ei and ED. As t — > ±00, 

|| *(., t, e) - i>°(;t, e,T(r)\\ LH ^) = O e (l/\tf), 

for any < j3 < 1/2. 

Remark 5.3. The estimate above is valid for any value of e > 0, and its proof 
essentially relies on integration by parts. 

We now analyze the freely propagating states ipj(x, t, e, ±). We do this by in- 
vestigating the small e behavior of the coefficients cJ(±oo, E, e). We are interested 
in solutions associated with only one electronic state in the remote past. For con- 
creteness, we choose level 2. Proposition 15.21 shows that this corresponds to the 
initial condition 

(5.7) cj(-cx),£,e) = S ji2 5 a< _. 

The component that has made the transition to level 1 in the course of the evolution, 
and propagates in the positive direction for large positive times is ipi (x, t, e, +). It 
is characterized by c 1 "(+oo, E, e). We analyze solutions of the equation satisfied by 
cj by means of complex WKB methods, along the lines of |24| . since we are in an 
avoided crossing regime. At this step, we require the parameter S > to be small, 
so that we can apply the results of |24| . They show that there exist So > 0, and 
8 e o(<5), such that for 5 < 5q and e < eo(S), 

(5.8) c^(+oo,E,e) = e -i»(C) e </ c MJS,*)<i*A 2 (i + B (e 2 )), for all E E A. 

Here £ is a loop in the complex plane, based at the origin, which encircles the zero 
z (S) E C \ R of k 1 (z,E) - k 2 {z,E) closest to the real axis. Also, 9(() E C only 
depends on the analytic continuation along £ of the electronic eigenvectors. 
We now have everything to state the main result of 21 : 

Theorem 5.4. Let ip{x,t,e) be a solution to the molecular Schrodinger equa- 
tion 14 -3\) with electronic Hamiltonian h(x, S) satisfying hypotheses M e ; and AC e i . 
Assume the solution is characterized asymptotically in the past for negative x's by 

lim |j ip(x,t,e) - ip2(x,t,e,-) \\wr-) =0, 
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with 

fc(x,t,e,-) = 4> 3 {X) I Q(^ £ ) e -itE/e* e i( X k 2 (-oo,E) + ^(E))/e* dE 

2 ' /a ^/2k 2 (-oo,E) 

where the energy density Q(E, e) satisfies ED. Let E(k) = k 2 /2 + i?i(+oo) 

a{E) = G(E) + Im J k 2 (z,E)dz, and 
C 

k{E) = J{E) -Ref k 2 (z,E)dz - uj+(E). 
C 



Assume E* is the unique absolute minimum of a(-) and set k* = ^2{E* - £1(00)). 

Then, there exist So > 0, p > arbitrarily close to 5/2, and a function eo : 
(0, So) — > M + , such that for all < [3 < 1/2, 5 < So, and e < £o(S), the following 
asymptotics hold as t — > oo: 

1pl(x,t,E,+) = 

-i0(C) e 3/2 3/4 -a(E*)/e 2 

^iW Tu— e 4S+(t)/e MA + (t),B + ,e",a + (t),n + ,x) 

(£a(E(k))\ k *) 1 

x P(E*,e) ^ e -^*)-fc* 2 «'(B*))/ £ 2 + 0( e -«(B*)A 2 e P) + O e (l/t"), 
where ipo is a Gaussian j5.4\ ) parametrized by 



V+ = k*, a+(t) = k'{E*) + k% B+ = l/^d 2 a(E{k))\ k=k , , 
A+(t) = B+ ( d 2 k a(E(k))\ k=k * + i(d 2 k K(E(k))\ k=k , + 1) ) , and 
= (fc* 2 /2-^i(+oo)) t. 

All error terms are estimated in the L 2 (R) norm, and the estimate 0(e _Q(ir )/e2 e p ) 
is uniform in t. The O t {l/t^) may depend on e. 

Remark 5.5. The leading order of the piece transmitted to level 1 is always 
given by an exponentially small prefactor times a freely propagating Gaussian nu- 
clear wave packet, i.e. a Gaussian semiclassical wave packet in the constant poten- 
tial given by the asymptotic electronic level. This is true even if the incoming state 
is not Gaussian. 

All quantities computed from the electronic Hamiltonian depend on 5, despite 
this not being emphasized in the notation. Although we shall not do it here, one 
can further investigate the characteristics of the asymptotic state for small S's. 
Without being too precise, we mention that if k c (E) — y/2(E — E 2 (0, 0)) is the 
classical momentum at the crossing point of the electronic levels for S = 0, one has 

a(E) = g(E - E ) 2 /2 + + 0((E - E ) 3 ) + O(o 3 ). 

The factor <5 2 7r/(4 k c (E)) is the Landau-Zener decay rate computed from the elec- 
tronic levels divided by the classical momentum at the avoided-crossing, and g > 
stems from the energy density. So, for small values of S and for a sufficiently nar- 
row energy window A, a is a quadratic term plus a positive, decreasing, convex 
function of E. This implies that E*, determined by a'(E*) = 0, is strictly larger 
than Eo and depends explicitly on the energy density. Furthermore, the decay rate 
has a{E*) < a(E ). This has the following consequences. 
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Remark 5.6. The average momentum k* is larger than what a naive applica- 
tion of energy conservation yields, and the exponentially small prefactor is larger 
than the prediction given by the Landau-Zener formula for h(a(t)). The correct 
values of both of these quantities depend explicitly on the chosen energy density. 

The lesson we learn from this is that the leading behavior of the piece of the 
wave function that describes transitions from one electronic level to the other cannot 
be determined by straightforward semiclassical considerations and the Landau- 
Zener formula. The way the incoming asymptotic state is prepared in the remote 
past plays a crucial role in the main characteristics of the asymptotics. 

Results of this sort are obtained in |26j for scattering systems defined by more 
general semiclassical autonomous PDEs entertaining dispersive waves. While the 
overall general strategy remains the same, the separation of variables, study of 
the stationary equation by means of complex WKB methods and analysis of the 
space-time properties of exact solutions, can be quite different. 

The description of the time development of the transitions between the elec- 
tronic levels represents a challenging question. Betz and Teufel have heuristic 
arguments and numerical evidence supporting a description of this phenomenon. 
However, the mathematical problem remains open. 
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